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Abstract. Motivated by examples of symmetrically constrained compositions, super con- 
vex partitions, and super convex compositions, we initiate the study of partitions and com- 
positions constrained by graph Laplacian minors. We provide a complete description of the 
multivariate generating functions for such compositions in the case of trees. We answer a 



O , question due to Corteel, Savage, and Wilf regarding super convex compositions, which we 

describe as compositions constrained by Laplacian minors for cycles; we extend this solution 



to the study of compositions constrained by Laplacian minors of leafed cycles. Connections 
are established and conjectured between compositions constrained by Laplacian minors of 
leafed cycles of prime length and algebraic/combinatorial properties of reflexive simplices. 



1. Introduction 



. , A partition of a positive integer n is a sequence A := (Ai, A2, . . . , Afc) of non-negative inte- 

. gers satisfying A^ < Afc_i < ■ ■ ■ < Ai such that J2i -^j = ^- We call each Aj a part of A and 

O ■ say that A has k parts. A (weak) composition of n is a sequence A := (Ai, A2, . . . , A^) of non- 

^ I negative integers satisfying Aj = n. Beginning with [Ij and through eleven subsequent 

O I papers, George Andrews and various coauthors successfully revived Percy MacMahon's tech- 

CN ■ 

nique of Partition Analysis (PA) from obscurity as a tool for computing generating functions 

for partitions and compositions. PA is particularly well-suited to the study of partitions and 

^ ■ compositions constrained by a linear system of inequalities, i.e., of the form AX > b for some 

H ! 

5^ , integral matrix A and integral vector b; this is equivalent to the study of integer points in a 

rational polyhedral cone. Motivated in part by this renewed interest in PA, there has been 
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growing interest in the application of polyhedral-geometric techniques to the study of integer 
partitions and compositions. 

For example, in [21], Igor Pak considers several partition theory problems from the point 
of view of lattice-point enumeration, leading to the establishment of unique bijective corre- 
spondences. In [1], Matthias Beck, Ira Gessel, Sunyoung Lee, and Carla Savage consider the 
general case of determining integer-point transforms of symmetrically constrained families 
of compositions; their techniques include both unimodular decompositions of cones and the 
use of permutation statistics. In [llj, Katie Bright and Carla Savage produce a bijection 
between certain "boxes of width t" of lecture hall partitions and integer points in the cube 
[0,^]"^; the methods for defining this bijection are motivated via geometry and again involve 
permutation statistics. In [5], the first author, Matthias Beck, and Nguyen Le give geometric 
proofs of several of the results in the PA papers by George Andrews et al; these proofs often 
produce multivariate generating functions in cases where PA produced only univariate iden- 
tities. Finally, motivated by [4j and [11], the first author and Matthias Beck use geometric 
techniques to establish new multivariate Euler-Mahonian distributions in [TU] . 

This paper continues this line of investigation by initiating the explicit study of partitions 
and compositions constrained by graph Laplacian minors. Algebraic graph theory is a rich 
source of interesting constraint matrices, the first being the signed incidence matrix of a 
graph, defined as follows. Let G = ({0, 1, . . . , n}, E) be a simple graph with n + 1 vertices 
and \E\= m edges. 

Definition 1.1. For any fixed edge orientation e, the signed vertex-edge incidence matrix 
of G with respect to e is the + 1 x m matrix d, with rows indexed by vertices of G and 
columns indexed by edges of G, whose {v,e)-th entry is 1 if v is the positive end of e, —1 if 
V is the negative end of e, and otherwise. 

The study of compositions with constraint matrices given by d'^ for various graphs was 
recently advanced in [13] by using five guidelines from MacMahon's partition analysis along 
with the Omega software package. Additionally, such compositions often fall within the realm 
of the theory of P-partitions [2H Chapter 4]. Our goal is to extend this line of investigation 
by using as constraint matrices various minors of graph Laplacians, defined as follows. 

Definition 1.2. The vertex Laplacian of G is L := D — A, where D is the n + 1 x n + 1 

matrix indexed by the vertices of G such that the matrix entry dij is the degree of vertex i 
if i = j (^iT-d otherwise and A is the 0/1-valued adjacency matrix of G indexed in the same 
way as D. Given a graph G, the i-th Laplacian minor of G is Li, the Laplacian matrix L 
with row i and column i deleted. 
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The connection between the incidence matrix of a graph and its Laplacian is the well- 
known equation L = dd^ [51 Chapter 6]. Kirchoff's Matrix- Tree Theorem [S] asserts that 
given a graph G on {0, . . . , n} and given any i between and n, det(Lj) is equal to the number 
of labeled spanning trees of G. This establishes a strong connection between Laplacian minors 
and enumerative properties of graphs. 

There are several examples of partitions and compositions occurring in the literature on 
partition and composition enumeration that are constrained by graph Laplacian minors. The 
following three such examples provide further motivation for the study of these objects; note 
that these examples were not originally introduced in the context of graph Laplacian minors. 

Example 1.3. Let Kn+i be the complete graph on {0, 1, . . . , n}. The compositions determined 
by LqA > are one example of the class of symmetrically constrained compositions. A study 
of multivariate generating functions for symmetrically constrained compositions was given 
by Beck, Gessel, Lee, and Savage in [1]. 

Example 1.4. Let Pn be the path graph with vertices {0, l,...,n — 1}. The partitions 
determined by LqX > are named super convex partitions by Corteel, Savage, and Wilf 
in Example 3 of [12]; they were first studied using techniques from MacMahon's Partition 
Analysis by Andrews in pp, where they are named partitions with mixed difference conditions. 
The multivariate generating function for such partitions is described in [T2] . 

Example 1.5. Let C„ be the cycle on vertices {0, 1, . . . ,n — l}. The compositions determined 
by LqX > were introduced by Corteel, Savage, and Wilf in the concluding remarks of [12] . 
where they are named super convex compositions. The determination of the multivariate 
generating function for super convex compositions was given as an open problem at the end 
of [ig. 

We now introduce our primary object of study. 

Definition 1.6. Given a graph G with Laplacian L and a vertex i of G, the compositions 
constrained by the i-th Laplacian minor of G are those satisfying LiX > 0. 

Our paper is structured as follows. In Section [21 we review techniques from integer point 
enumeration in polyhedral cones that will be useful for producing generating function iden- 
tities. In Section [21 we investigate compositions constrained by Laplacian minors for trees. 
Theorem 13.11 gives a combinatorial interpretation for the entries of the matrix inverse for 
tree Laplacian minors, leading to Corollary 13. 2^ the main result of the section. 

In Section [H we provide a solution to the open problem of Corteel, Savage, and Wilf 
discussed in Example 11.51 However, as shown in the section, our solution is less clear than 

3 



one might want. Thus, motivated by the underlying graphical structure of our compositions, 
in Section Owe examine compositions constrained by Laplacian minors for leafed cycles, i.e. 
for cycles with an additional vertex forming a leaf appended. These compositions are closely 
related to the super convex compositions in Example II. 5[ The addition of a leaf to the cycle 
causes a change in the structure of the resulting compositions, leading to Theorem 15.11 We 
also consider in Conjecture 15.51 a possible combinatorial interpretation for the compositions 
constrained by leafed cycles, and we prove special cases of this conjecture. 

In Section |6l we consider generating functions for certain subfamilies of compositions con- 
strained by Laplacian minors for leafed cycles of prime length. Surprisingly, these generating 
functions are Ehrhart series for a family of reflexive simplices, as shown in Theorem 16.31 Af- 
ter introducing Ehrhart series and reflexive polytopes and proving Theorem 16.31 we indicate 
some possible connections between our enumerative problem and the algebraic/combinatorial 
structure of these simplices. 

We close in Section [7] with a conjecture regarding compositions constrained by Laplacian 
minors for leafed cycles of length 2^. 

2. Cones and generating functions 

The problem of enumerating integer compositions and partitions satisfying a system of 
linear constraints AX > b is identical to that of enumerating integer points in the polytope 
or polyhedral cone defined by these constraints. For further background related to integer 
point enumeration in polyhedra, see [B]. 

Given a square matrix A with det(74) 7^ 0, the solution set 

C = {x eW : Ax >0} 

is called an n- dimensional simplicial cone. The reason for the terminology "simplicial" 
comes from C being an ra- dimensional cone with n bounding hyperplanes, also called facets. 
The following lemma, whose proof is straightforward, indicates that simplicial cones can be 
described either by specifying these n facets or by specifying the n ray generators for the 
cone. 

Lemma 2.1. Let A be an invertible matrix in M"^". Then 

{xeW : Ax>0} = {A-^y : y > 0} 

where each inequality is understood componentwise. 
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We encode the integer points in a cone C using the integer point transform of C given by 
the multivariate generating function 



m2 _ _ ^m„ 



This function encodes the integer points in C as a formal sum of monomials. Throughout 
this paper, we will typically need to specialize ac in one of two ways. If we set zi = Z2 = 
■ ■ ■ = Zn = q, we obtain 



where \m\ = mi + 1712 + • — h m„. Setting zi = q and Z2 = z^ — ■ ■ ■ — Zn — 1, we obtain 



aciq,l,...,l)^ac{q)^ J2 



meCnZ" 

When A is an integer matrix having det{A) — 1, the ray generating matrix for the n- 
dimensional simplicial cone C = {x e : Ax > 0} has integer entries; hence, every integer 
point in C can be expressed as A~^y for a unique integer vector y > 0. In this case, the 
integer point transform can be simplified as a product of geometric series as follows. 

Lemma 2.2. Let A — {aij) be an integer matrix with det(A) = 1 and inverse matrix 
A-^ = (bij). Then 

1 



ae{zi,Z2, ...,Zn):-- 



More generally, suppose det(A) 7^ 1. In this case, ac{zi, Z2, ■ ■ ■ , Zn) is still rational, but 
need not have integer entries. As a result, representations of integer points in C are not 
necessarily obtained by using integral combinations of the minimal integral vectors on the 
generating rays of C. We can work around this difficulty in the following manner. 

Let f 1, V2, ■ ■ . ,Vn be the (integral) ray generators of C. Define the (half-open) fundamental 
parallelepiped (FPP) to be 

n = {Xivi + X2V2 H XnVn : < Ai, A2, . . . , A„ < 1} . 

The cone C can be tiled by integer translates of 11, but there may be integer points contained 
in n that can only be obtained using strictly rational scaling factors of the ray generators of 
C. This leads to our main tool for computing generating functions. 
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Lemma 2.3. Using the preceding notation, if the points contained in H are wi, W2, ■ ■ ■ , Wm, 
the integer point transform of C is given by 

i=l ^1 ^2 ■ ■ ■ 



(1) crc{zi, Z2, . . . , Zn) := 



'■J 

3. Trees 

We begin our study of compositions constrained by graph Laplacian minors with the 
special case of trees. Since a tree T has a unique spanning tree, namely itself, the Matrix- 
Tree Theorem implies that det(Lj) = 1 for any vertex i. Thus, compositions constrained by 
tree Laplacian minors will have a generating function of the form 

1 



ac{zi,Z2, ...,Zn) 



Since the &fcj's arise as the entries in the inverse matrix for Li, we are interested in calculating 
this inverse. We first observe that without loss of generality we may consider only those 
vertices i which are leaves of T, through the following reduction. 

Suppose that a vertex i G T is not a leaf. Then taking the minor at i deletes row and 
column i from the Laplacian matrix, which equates to removing vertex i in the graph. We 
note that removing the vertex i will split our tree into some number of disjoint subtrees, 
which we label Si, . . . , Sk- Label the vertices of Si consecutively from ri + 1, ri + 2, . . . , r2, 
then label the vertices of S2 consecutively as r2 + 1, r2 + 2 . . . , rs, and so on until labeling the 
vertices of as + 1, . . . , r^+i. Now consider the r^+i x r^+i matrix M, indexed naturally. 
M has a block form resulting from the fact that the subtrees Si,...,Sk are distinct and 
disconnected. In particular, each subtree Sj, 1 < j < k, will have a corresponding block 
matrix Mj consisting of rows rj + l,..., r^+i and columns + 1, . . . , rj+i. 



M = 



Ml 







Mo 







Further, each Mj will correspond to the Laplacian minor for sj where the minor is taken 
at i, which is now simply a leaf of Sj. By properties of block matrices, we have that will 
be a block matrix with blocks . Thus, to find M~^, we compute the Laplacian minor 
inverse for each subtree. With this observation, we can begin our search for a general form 
of the generating function of compositions constrained by Laplacian minors of trees. 

3.1. Inverses of Laplacian minors for trees. To determine the denominator of ac for 
an arbitrary tree T minored at a leaf, we first find a combinatorial interpretation for the 
columns of the inverse of the Laplacian minor for T. (While Theorem 13. II below seems likely 
to be known to experts in algebraic graph theory, we could not find an explicit statement of 
it in the literature.) Throughout this section, we assume that T is a tree on the vertex set 
{0, . . . , n} with an arbitrary orientation of the edges of T. 

Theorem 3.1. Consider a tree T and its corresponding Laplacian matrix L. Minor at a 
leaf n ofT and let L^^ = {l-ij). Then entry lij is the distance from n to the path connecting 
the vertices i and j . Equivalently, lij is the length of the intersection of the path from n to 
i with the path from n to j. 

Theorem 13.11 has the following implication on the level of generating functions. Let ac{q) 
denote the specialiazed generating function ac{q, q, . . . ,q). 

Corollary 3.2. The specialized generating function for a tree T minored at leaf n is given 
by ^ 

where hi is the sum of the distances between n and the paths connecting vertex i with the 
other vertices of T . 

To prove Theorem 13. H we return to the fact that L = dd^ and proceed combinatorially. 
We say the n-th subminor of d (resp. of 9^) imposed by a vertex of a tree is the square 
matrix in which the n-th row, (resp. column) is deleted. We denote the n-th subminor of 
d as dn, with the understanding that (9j is (dn)'^. It is easy to see that since L = dd^ , we 
have Ln = <9„9^; thus, = (9^)"^9^^ 

Proposition 3.3. Let Gn = (gej), < j < n — 1, where Qej is determined as follows: 

• I if edge e is on the path between j and n and e points away from n; 

• —1 if edge e is on the path between j and n and e points toward n; 

• if edge e is not on the path between j and n. 
Then Gndn = I, hence (9^^ = Gn- 
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Proof. Set Gndn = M = {rriij) and first consider mj^j. Entries of row i of dn are non-zero for 
each edge, i.e., column, adjacent to i. The path connecting vertices i and n overlaps with 
exactly one of the edges adjacent to i, which is reflected in Gn- By construction, the entries 
in dn and G„ corresponding to this vertex/edge pairing have the same sign. So mj j = 1 for 
all i. 

Now consider rriij with i ^ j. The path from j to n overlaps with either no edges or 
exactly two edges connected to i. If there is no overlap, then the product of row i in and 
column j in dn is zero. If there is overlap, suppose both edges point either toward or away 
from n. Then the entries in G„ have the same sign, but those in dn have different signs, thus 
summing to zero. A similar canceling occurs if the edges point in different directions. So, 
rrii^j = if i 7^ j, implying M = I. Therefore, On = d~^. □ 

Proof of Theorem \3.1[ Since = {dn)~^dn^ , we may multiply columns i and j in 
component- wise to find /j.-,-. Consider an arbitrary edge of the graph. If the edge belongs to 
only one of the paths from i to n or from j to n, then there will be a zero corresponding to 
that edge in either column i or j, contributing nothing to entry lij. Otherwise we obtain 
a 1, since the sign of a row for any edge is constant. Thus, summing the component- wise 
product will record the length of the path from n to the path between i and j. □ 

Theorem 13.11 provides some insight into the form of the generating function for composi- 
tions constrained by an arbitrary tree. However, determining the lengths of paths between 
vertices in trees is not an easy exercise in general. In the next subsection, we analyze a 
special case. 



3.2. k-ary trees. We next provide an explicit formula for a specialized generating function 
for compositions constrained by Laplacian minors of k-aij trees. Recall that a fc-ary tree 
is a tree formed from a root node by adding k leaves to the root, then adding k leaves to 
each of these vertices, etc., stopping after a finite number of iterations of this process. To 
study the compositions constrained by a Laplacian minor of a A;-ary tree T, we add a leaf 
labeled to the root of T and minor at this new vertex. We define the level of a vertex 
as its distance to vertex and a subtree of a vertex as the tree which includes the vertex 
and all of its children (we assume all edges are oriented away from vertex 0). Recall that 
[n]g ■=l + q + q^ + q^ + --- + 

8 



Theorem 3.4. For a k-ary tree with n levels, 

1 



JJ A _ qj[n - 3 + + k-^m]k - 1) - - ^))\ 

j=i V / 



Proof. Suppose a k-ary tree T is given. In order to find the generating function consider the 
column sums of Lq ^. Choose a particular column which will correspond to a vertex v E T. 
The column sum will be the sum of the distances between and all distinct paths in the 
tree with v as an endpoint. By the symmetry of k-ary trees, we choose an arbitrary vertex 
V at level j e [n] , since the corresponding column sum will be the same as any other vertex 
on that level. 

We claim that j[n — j + l]k + '^IZiU ~ i){k'^~^^~'^^) is the column sum corresponding to 
any vertex at level j. For any vertex c which is a child of v (and including v), the distance 
between and the path from v to c will be j. There are 

1 + A; + A;' + • • • + A;""^' = = [n - j + l]k 

K 1 

such vertices and each of these vertices contributes j to the column sum. Hence in total these 
vertices contribute the first term to our claimed exponent. 

Consider the path between and P = {v = Wo, Wi, W2, . . . , Wj-i, 0}. For each i G [j — 1], 
define Si to be the set of vertices which are in the subtree with Wi as the root but not in the 
subtree with Wi-i as the root. The vertex Wi is on level j —i of T. Thus, the subtree whose 
root is Wi has n — {j — i) levels of its own. So, there are + + ■ ■ ■ + k^~^^~'^^ vertices in 
the subtree with Wi as a root, and A;" + A;^ + • • • + k^~^^~^^~^ vertices in the subtree of Wi-i. 
Thus, there are 



A;° + A;^ + • • • + A;"-^^-^) - (A;° + A;^ + • • • + A;"-^^'-^)-^) = A;"" 



vertices in Si. By definition, the shortest distance between and the path connecting v and 
any vertex s E Si is j — i, since the path connecting v and s passes through Wi, which must 
be the closest point on this path to 0. 

Therefore, for each i E [j — 1], the sum of the columns corresponding to vertices in Si will 
be (j — i)(A;"~*^-'~*)). So, for all vertices which are not children of v, the total column sums 
will be jyiZlij — i){k^~^^~^^). Hence, j[n—j + l]k + '^{ZliJ ~ i){k'"~^^~^^) is the column sum 



corresponding to any vertex at level j G [n]. Since there are y ^ vertices on each level, for 
a fixed j, the term 

1 _ qj[n-J + l]k + k--'iji[j]k - 1) - kiilj], - 1)) 

will appear times as an exponent in the denominator of our generating function. □ 
When k = 2, the formulas for our exponents simplify nicely, as follows. 



Corollary 3.5. For a binary tree with n levels, 
crciq) = 



1 



J](l-gJ(2"-^■-^^-l) + E^:^(J-0(2 
j=i \ 

1 



23 



n(l-g2"-^H2^-l)-j) 



2^' 



-1 



We omit the proof of the equivalence of the two rational forms in our Corollary, as they 
are the result of a tedious calculation. For those interested in the calculation, we mention 
that the key step is the use of Et!(j - ^2' = Eti^(2 + 2^ + ■ ■ ■ + 2^). 



4. ra-CYCLES 

In this section we consider the open problem due to Corteel, Savage, and Wilf given in 
Example II. 5[ Recall that the Laplacian minor for a general n-cycle on vertices 1, 2, . . . , n, 
labelled cyclically and minored at n, has the following form when rows and columns are 
indexed 1, 2, . . . , n — 1. 



/ 2 -1 
-12-10 
0-12-1 



-'n,cyc 





V 


















Define 



\ 





-1 2 -1 
0-12/ 



{A G R" ■■ Ln,cyc>^ > 0} and (T^„(zi, . . . , z„_i) := 
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AeK:„nZ"-i 



Our goal is to study the integer point transform for /C„ by enumerating the integer points 
in the fundamental parallelepiped for /C„; along the way, we will see that there is a hidden 
additive structure to the ray generators of )Cn which is useful. 

By Lemma [2. H the rows of Ln^cyc determine the normal vectors to the facets of /C„, while 
the columns of L~^^y^ provide (non- integral) ray generators of /C„. We first explicitly describe 
the entries in L~^^y^. 

Proposition 4.1. For an n-cycle and 1 < i, j < n — 1, L~\ = (bij) where 



i(n—j) t ■ ^ ■ 
¥ J 

j(n—i) -e ■ 

— if l> ] 



n,cyc 



Proof. Define the matrix B = (bi j) using the 6- values given above. Let B ■ L 
since Ln,cyc and B are symmetric matrices, it follows that ajj = aj^i. Entry aij is given by 
26j J — bi^x — bi^y where x and y are the neighbors of j in the cycle, where we drop the 6-terms 
for X or y equal to n. It is straightforward to show that BLn^cyc = I through a case-by-case 
analysis. 

As a representative example, for i < j we have 



26,j- — bij-i — bn,j+i 



i[n- ]] 



n 



i{n- j + 1) 



n 



i{n- j - 1) 



n 



□ 



We next scale the matrix Lj^].y^ by n so that the columns are integral; in doing so, we find 
that there is group-theoretic structure to the ray generators of /C„, in the following sense. 

Proposition 4.2. Let the matrix M := nL~\ (mod n). Then 



M 



Vi V2 ■■■ Vn-2 Vn-1 



where f i = (n — 1, n — 2, . . . , 2, 1)^ G ZJJ ^ and Vk = kvi e 



n-l 
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Proof. Let M = {rriij). Since rriij is equal to either i{n — j) or j{n — i), it follows that 
rriij = —ij mod n. Consider column Vi. We find 

rp 

Vi = (—1 ■ 1, — 2 ■ 1, . . . , —{n — 1) ■ 1) mod n, 
= {n-l,n-2,...,lf eZ^^~\ 



Evaluating an arbitrary column Vk, 1 < k < n — 1, we have 

Vk = (—1 ■ k, —2 ■ k, . . . , —{n — 1) ■ kf" mod n, 

= k-vie zi-\ 

□ 

Our next step is to use the columns of nL~^^y^ as ray generators for /C„, leaving us only to 
understand the integer points in n„, the fundamental parallelepiped for /C„. For any integer 
point A G /C„, since L^^cyc is an integer matrix, L^^cyc^ = c G Z"^^ Thus, every integer point 
in /C„ can be expressed in the form L~^^y^c for an integral c; further, if A G n„, then 



where the entries in c are necessarily taken from {0, 1, 2, ... ,n — 1}. The group-theoretic 
structure of nL~^^y^ (mod n) allows us to classify these points. 

Proposition 4.3. Suppose that the point x = {xi,X2, ■ ■ ■ ,Xn-i) G n„ is of the form L~\y^c 
for some vector c with entries taken from {0, 1, 2, . . . , n - 1}. Then 
/C„ is an integer point if and only if xi is an integer. 

Proof. If X is an integer vector, then xi is certainly an integer. For the converse, suppose that 
Xi is an integer. This happens precisely in the case where, for I denoting the first row of nL~^, 
we have that l-c = X]j'=i an integer multiple of n. Thus, Vi-c = Yl]=i '^ij'^i = niod n 

where fi is as defined in Proposition 14. 2[ Since Vk = kvi, 

n n n 

VkjCj = kvijCj = k VijCj = mod n. 
j=i j=i j=i 

Therefore, Xk is also an integer. □ 
Thus, the integer points in are parametrized by the set of solutions to the system 

n-l 

— j)cj = (mod n) 
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where Cj G {0, 1, 2, . . . , n — 1}. Equivalently, these points are parametrized by the set of 
partitions of multiples of n into positive parts not exceeding n — 1, with no more than n — 1 
of each part. The solution set to this system of equations has arisen before in the study of 
"Hermite reciprocity" given in |T5j. Our next theorem follows directly from the preceding 
discussion. 

Theorem 4.4. Let Sn denote the set of solutions to the system X]j'=i(^ ~ = (modn) 
where Cj G {0, 1, 2, . . . , n — 1}. Let denote the i,j-th entry of nL~^^y^. Then 



While Theorem 14.41 and its specializations resolve the question of Corteel, Savage, and 
Wilf, this resolution is not as clean as one might like. In the next section, we show how a 
slight modification of the graph underlying this problem leads to much more simply stated 
and elegant, but still closely related, univariate generating function identities. 



In this section, we examine the case of compositions constrained by n-cycles augmented 
with a leaf; we will refer to such a graph as a leafed n-cycle. Throughout this section, we 
label the leaf vertex as n and label the vertices of the n-cycle cyclically as 0, 1, . . . , n — 1, 
where n is adjacent to as depicted in Figure [H 




5. Leafed n-CYCLES 



1 



2 




n-1 



n-2 



Figure 1. A leafed n-cycle 
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For a leafed n-cycle, we denote by L„ the minor of the Laplacian taken at vertex n and 
observe that it has the form 

/ 3 -1 ■■■ 
-12-10 ••• 
-1 2 -1 ••• 
-1 2 -1 ... 



V 






























-l\ 







-1 





where the columns (resp. rows) are labeled from top to bottom (resp. left to right) by 



0,1,2,. 



1. Define 
C 



{AgM":L„A>0} and cJcM ■= Yl ^' 

AGC„nZ" 

While not immediately apparent, the reason for focusing on the first coordinate of integer 
points in C„ will become clear. 

5.1. A generating function identity. Our main result regarding leafed cycles is the fol- 
lowing generating function identity. 

Theorem 5.1. Let Sn denote the set of solutions to the system 

n-l 

■ Co + "Y^in — j)cj = (mod n) , 



where Cj G {0, 1, 2, . . . , n — 1}. For c G Sn, define 



Then 



E 



Note that YlceSn ^'^^'^^ ordinary generating function for the "number of parts" statistic 

on the set of all non-negative partitions of multiples of n into parts not exceeding n — 1, with 
no more than n — 1 of each part. The similarity between Theorem 15.11 and Theorem 14.41 is 
immediate. However, by switching from ra-cycles to leafed n-cycles, a more elegant univariate 
identity is obtained than those obtained by straightforward specializations of Theorem 14.41 
We remark that our proof of Theorem 15.11 can easily be extended to produce a multivariate 
statement analogous to Theorem 14. 4( we omit this, since it is the specialization above that 
we find the most interesting consequence of the multivariate identity. 

As in the previous section, our goal is to apply Lemma 12.31 We first explicitly describe 
the entries in L~^, which are closely related to the entries of L~^^y^. 
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Proposition 5.2. For a leafed n-cycle and < i, j < n — 1, L^^ — (bij) where 



iin-j) 
n 



+ 1 Zfl<j 

+ 1 tfi>j 



"1,3)1 



: since 



Proof. Define the matrix B — {bij) using the 6- values given above. Let LnB - 
L„ and B arc symmetric matrices, it follows that For i 0, entry aij is given by 

26i 1 — bxj — byj where x and y are the neighbors of i in the cycle; the case i = is similar. 



From these observations, it is straightforward to show that LnB — I through a case-by-case 
analysis. 

For example, for j ^ 0, qqj = is shown as follows. 



^^0,7 — 3^0,1 ~ bii — b. 



3 
0. 



O(n-j) 



n-l,j 



+ 1 - 



n 



l{n-j) 



+ 1 



n 



j{n - (n - 1)) 



+ 1 



n 



□ 

Note that every entry in the top row (and in the first column, by symmetry) of is 
equal to 1; it is this property of the inverse that allows us to effectively study the generating 
function recording only the first entry of each composition in C„. We next scale the matrix 
by n so that the columns (i.e. ray generators of our cone) are integral; in doing so, we 
find that there is group-theoretic structure to the columns of L~^, again related closely to 
that for L-^^y^. 



Proposition 5.3. Let the matrix M — nL„^(mod n). Then 



M 



\ 



Vo Vi V2 



Vn-2 Vn-l 



J 



where i^i = (0, n — 1, n — 2, . . . , 2, 1)'^ e and Vk — kvi e Z^. 



Proof Let M 
that rriij 



). Since rriij is equal to either i{n — j) + n or j{n — i) + n, it follows 



-ij mod n. Consider column vi. We find 

vi = (-0 • 1, -1 • 1, -2 • 1, . . . , -{n - 1) • 1)^ mod n, 
^{0,n-l,n-2,...,lf eZ^. 
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Evaluating an arbitrary column Vk, < k < n — 1, we have 



Vk = (—0 ■ k, —1 ■ k, —2 ■ k, . . . , — (n — 1) ■ k)^ mod n, 



= k-vi eZ: 



□ 



As in the case of n-cycles, our next step is to use the columns of nL~^ as ray generators 
for Cn', we must understand the integer points in n„, which now denotes the fundamental 
parallelepiped for C„. For any integer point A G C„, since L„ is an integer matrix, L„A = c G 
Z". Thus, every integer point in Cn can be expressed in the form L~^c for an integral c. As 
before, if A G n„, then 



where the entries in c are necessarily taken from {0, 1,2, ... ,n — 1}. The group-theoretic 
structure of nL~^ (mod n) allows us to classify these points. 

Proposition 5.4. Suppose that the point x = {xq, Xi, . . . , x„„i) G n„ is of the form L~^c for 
some vector c with entries taken from {0, 1, 2, . . . , n — 1}. Then x = {xq, Xi, . . . , Xn-i) G C„ 
is an integer point if and only if Xi is an integer. 

Proof. The proof is almost identical to that for the n-cycle case. □ 
Thus, the integer points in n„ are parametrized by the set of solutions to the system 



where cj G {0,1,2, ... ,n — 1}. Equivalently, these points are parametrized by the set of 
partitions of multiples of n into non- negative parts not exceeding n — 1, with no more than 
n — 1 of each part. 

Our final observation in this subsection, and the key observation needed to apply Lemma [273] 
and complete the proof of Theorem 15.11 is that for any point A such that L„A = c G n„, we 
have that Aq = X^i^i- This is a consequence of Proposition 15.21 which shows that the first 
row (and column) of are vectors of all ones. 

5.2. Cyclically distinct compositions. This subsection is devoted to the following con- 
jecture, which provides a combinatorial interpretation for the coefficients of ac„{q). 

Conjecture 5.5. The number of integer points A G €„, with Xq = m is equal to the number 

of compositions of m into n parts that are cyclically distinct. 
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n-l 




Conjecture 15.51 is based on experimental evidence for small values of n; further, we used 
LattE [TS] to check the conjecture for all prime values of up to = 17. An obvious way 
to attack Conjecture 15.51 is to attempt to prove that for any A G C„ with Aq = m, we have 
that A = L~^c for an integer vector c whose entries form a weak composition of m that is 
cyclically distinct from all other such c vectors forming elements of C„. Unfortunately, this 
does not work out, as the following example demonstrates. 

Example 5.6. The weak compositions of 3 into 3 parts are given by the following families 
Fi,F2,F3, and F^; each family is the orbit of a single composition under the cyclic shift 
operation. 

Fi: (3,0,0) (0,3,0) (0,0,3) 
F2: (2,1,0) (0,2,1) (1,0,2) 
F3: (1,2,0) (0,1,2) (2,0,1) 
F4: (1,1,1) 

The four compositions which satisfy L^^c E Cn are the elements of Fi and F4. Thus, while 
this shows that Conjecture 15.51 holds in this case, it is not a consequence of the c- vectors 
being themselves cyclically distinct. 

While this obvious approach to Conjecture 15.51 fails in certain cases, it is successful in 
others, as seen in the following theorem. 

Theorem 5.7. Let n he prime and fix m such that gcd(m, n) = 1. The number of integer 
points in Cn with first coordinate m is equal to the number of compositions of m into n parts 
that are cyclically distinct. 

Proof. We must show that for every composition c of m with n parts, exactly one of the 
cyclic shifts of c will produce an integer point when is applied. Begin with a composition 
c = (co, Ci, . . . , Cn-i) of m such that L~^c G C„ fl Z". Then for all j, 

n-1 

CiVj^i = (mod n) 

1=0 

where Vj = {vj^, Vj^i, . . . , f j>-i) is the jth row of nL^^ ( mod n). Now consider another cyclic 
ordering c' of c where we shift the position of each entry by k, so = (cfc, Ck+i, ■ ■ ■ , Ck+n-i) 
with the indices taken mod n. Then the jth coordinate of L~^c' is given by 

n—l n—1 

i=0 i=0 
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By properties of vj given in Proposition |531 shifting the ith coordinate of Vj by k is equivalent 
to subtracting the row Vk from the row Vj component-wise. So, 

n— 1 n— 1 

1=0 i=0 

n— 1 n— 1 

ra-1 

= - fc ^ Q 

Since X^j'^To^ Q ^ mod n by assumption, and since Vj^k ^^=o Q ^ mod n for j,k ^ due 
to n being prime, we have that 

n-1 

Cit;^- ^ . 

Thus, the non-trivial cyclic shifts of c do not correspond to integer points in C^. □ 

It is straightforward to show that for the case where n is prime and n\m, for any composi- 
tion c of m into n parts that satisfies L~^c G C„ fl Z" all the cyclic shifts of c also satisfy this 
condition. Thus, the remaining challenge regarding Conjecture 15.51 in the case where n is 
prime is to show that this union of orbits of compositions under the cycle action contains as 
many compositions as there are orbits of the cycle action. While we have not been successful 
in proving this, the study of these compositions reveals a beautiful underlying geometric 
structure for our cone, which we discuss in the following section. 

6. Reflexive polytopes and leafed cycles of prime length 

For leafed cycles of prime length, the remaining cases of Conjecture 15.51 turn out to be 
equivalent to the problem of enumerating points in integer dilates of an integer polytope, 
which is the subject of Ehrhart theory. The polytope in question is a simplex, and we will 
prove that it has the additional structure of being reflexive, a property we now introduce. 

6.1. Reflexive polytopes. An n-dimensional polytope P in M" may equivalently [251 Chap- 
ter 1] be described as the convex hull of at least n + 1 points in or as the intersection 
of at least n + 1 linear halfspaces in M". When P = conv(V) and V is the minimal set 
such that this holds, then we say the points in V are the vertices of P. If P = conv(y) 
and |\^| = n -|- 1, then we say P is a simplex. We say P is integral if the vertices of P are 
contained in Z". 
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For P C an integral polytope of dimension n, and for t G Z>o, set tP := {tp : p G P} 
and Lp{t) := |Z" fl tP|, i.e. the number of integer points in tP is Lp{t). Define the Ehrhart 
series for P to be Ehrp(a;) := 1 + 'Ylit>i Lp{t)x* . A fundamental theorem due to E. Ehrhart 
states that for an n-dimensional integral polytope P in M", there exist complex values 



Ehrp(x) = and \" h* ^ 



X] 

3 



h* SO that 



A stronger result, originally due to R. Stanley in |23], is that the hy?, are actually non- 
negative integers; Stanley's proof of this used commutative algebra, though several com- 
binatorial and geometric proofs have since appeared. We call the coefficient vector h*p := 
(/iq, . . . , h*^ in the numerator of the rational generating function for Ehrp(a;) the h-star vector 
for P. The volume of P can be recovered as (Ylij h*)/n\. 

Obtaining a general understanding of the structure of /i*-vectors of integral polytopes is 
currently of great interest. Recent activity, e.g. [21 [5l [TJ HI [T71 [161 1201 122], has focused on 
the class of reflexive polytopes, which we now introduce as they will be needed in Section [61 
Given an n-dimensional polytope P, the polar or dual polytope to P is 

= {x G M" : a; ■ p < 1 for all p G P} . 

Let P° denote the topological interior of P. 

Definition 6.1. An n-dimensional polytope P is reflexive ifOEP° and both P and P^ are 
integral. 

Reflexive polytopes have many rich properties, as seen in the following lemma. 

Lemma 6.2. [21 [IT] P is reflexive if and only if P is an integer polytope with G P° that 
satisfies one of the following (equivalent) conditions: 

(1) P^ is integral. 

(2) Lpo(t + 1) = Lp{t) for all t G Z>o, i.e. all lattice points in sit on the boundary of 
some non-negative integral dilate of P. 

(3) h* = for all i, where h* is the i*'^ coefficient in the numerator of the Ehrhart 
series for P. 

Reflexive polytopes are simultaneously a very large class of integral polytopes and a very 
small one, in the following sense. Due to a theorem of Lagarias and Ziegler [T^], there are 
only finitely many refiexive polytopes (up to unimodular equivalence) in each dimension. 
On the other hand, Haase and Melnikov j^H] proved that every integral polytope is a face 
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of some reflexive polytope. Because of this "large yet small" tension, combined with their 
importance in other areas of mathematics as noted before, whenever a reflexive polytope 
arises unexpectedly it is a surprise, an indication that something interesting is happening. 

6.2. Reflexive slices of Laplacian minor cones. Our main result in this subsection is 
Theorem 16. 3[ which asserts that for prime n values, reflexive polytopes arise as slices of the 
cone constrained by Laplacian minors for leafed n-cycles. Note that when n is not prime, this 
construction does not always yield reflexive simplices; this has been verified experimentally 
with LattE [18] and Lemma 16.21 However, exactly when the construction described below 
produces reflexive polytopes is at this time not clear. 

Let p be an odd prime and let ICPp be the simplex obtained by intersecting the cone 
constrained by a leafed p-cycle with the hyperplane Ai = p in W. 

Theorem 6.3. ICPp is reflexive (after translation by an integral vector). 

Proof. By elementary results about polar polytopes found in [25, Chapter 2], P is a reflexive 
polytope if and only if P is integral, contains the origin in its interior, and has a half-space 
description of the form 

{xeR'^ : Ax> -1} 

where A is an integral matrix and 1 denotes the vector of all ones. 

We first observe that ICPp is clearly an integer polytope, as its vertices are the columns 
of pL~^, which are integral. A halfspace description of ICPp is given by Aq = p and 
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-1 


■ ■ 


■ -1 
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-1 ■ ■ 
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Ai 
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2 ■ ■ 


■ 




A2 
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■ ■ 


■ 2 




Ap-i 








Let Wi be the ith column of L^^, and observe that Lp(Yl,^Wi) = Lp(Yl,^L~^{ei)) = 1. Thus, 
Y2 - Wi is an interior integer point of our simplex ICPp] we want to change coordinates so that 
Y2- Wi is translated to the origin. We therefore consider solutions to our halfspace description 
of the form A + w,. It follows that the description of our translated ICPp does not use a 
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Ao-coordinate (since Xq + p = p implies Aq = 0), and is given by 
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Since ICPp is an integral simplex, and Yli is an integer vector by which we translated ICPp, 
our translated ICPp is still integral. Thus ICPp is reflexive, and our proof is complete. □ 

While the proof above is elementary, it ties the study of compositions constrained by graph 
Laplacian minors into an interesting circle of questions regarding reflexive polytopes. The 
only remaining case of Conjecture 15.51 is that the number of integer points in m ■ ICPp is 
equal to the number of compositions of mp with p parts, up to cyclic equivalence, which is 
now asking for a combinatorial interpretation for the Ehrhart series of a reflexive simplex. 
Because of reflexivity, the generating function for this series yields a rational function with a 
symmetric numerator, which yields a functional relation on the Ehrhart polynomial for this 
simplex. 

The study of reflexive polytopes goes hand in hand with the study of several other in- 
teresting classes of polytopes; one such example are normal polytopes, where an integer 
polytope P is normal if every integer point in the m-th dilate of P is a sum of exactly m 
integer points in P. In our attempts to prove the prime case of Conjecture 15. 5[ we noticed 
that our techniques (though unsuccessful at providing a proof) provided evidence suggesting 
that ICPp is normal. Normality is implied by the presence of a unimodular triangulation for 
an integral polytope; while we are not yet certain of the existence of such a triangulation 
for ICPp or of the normality of ICPp, it would not surprise us if such a triangulation can be 
found for all prime p. 

Our final remark regarding ICPp, reflexivity, and normality, regards the unimodality of the 
h* vector of ICPp. It is a major open question (see [22] and the references therein) whether 
all normal reflexive polytopes have unimodal /i*-vectors. For p < 7, we found using LattE 
^18j that the /i*-vector for ICPp is unimodal. We again suspect that this holds in general. 

7. Leafed cycles of length a power of two 

We conclude our paper by considering leafed cycles for non-prime values of n; specifically, 
we study when n = 2^^ for some k. Using LattE [18], one observes that the cone constrained 
by a leafed 8-cycle does not have a reflexive 8-th slice (we suspect this to be the case for all 
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non- prime values of n). Nevertheless, the integer point transform of this cone exhibits some 
interesting "near-symmetry" in the numerator for small powers of 2. Based on experimental 
evidence, we offer the following conjecture. 

Conjecture 7.1. Let be the cone constrained by the Laplacian minor of a leafed n-cycle 
where n = 2^ for some integer k >2. Then the generating function has the form 




where f{q) has the following property. Let (ao, • • • ^dj) denote the coefficient list of f{q). If 
we append aQ to the end of this coefficient list, then take the difference between the appended 
coefficient list and its reverse, we obtain the coefficient list of the polynomial 

f ("-2)(_1)V» 

The near-symmetry of the numerator polynomial indicates that there should be some 
interesting structure to the finite group obtained by quotienting the semigroup of integer 
points in by the semigroup generated by a specific choice of integral ray generators for 
the cone. 

References 

[1] George E. Andrews. MacMahon's partition analysis. I. The lecture hall partition theorem. In Mathe- 
matical essays in honor of Gian-Carlo Rota (Cambridge, MA, 1996), volume 161 of Progr. Math., pages 
1-22. Birkhauser Boston, Boston, MA, 1998. 

[2] Victor V. Batyrev. Dual polyhedra and mirror symmetry for Calabi-Yau hypersurfaces in toric varieties. 
J. Algebraic Geom., 3(3):493-535, 1994. 

[3] Matthias Beck, Benjamin Braun, and Nguyen Le. Mahonian partition identities via polyhedral geometry. 
to appear in Developments in Mathematics, (memorial volume for Leon Ehrenpreis, edited by H. Farkas, 
R. Gunning, M. Knopp, and B. A. Taylor). Preprint at http://arxiv.org/abs/1103.1070. 

[4] Matthias Beck, Ira Gcsscl, Sunyoung Lee, and Carla Savage. Symmetrically constrained compositions. 
Ramanujan J.. 23:355 369, 2010. 

[5] Matthias Beck and Serkan Ho^ten. Cyclotomic polytopes and growth series of cyclotomic lattices. Math. 
Res. Let., 13(4):607-622, 2006. 

[6] Matthias Beck and Sinai Robins. Computing the continuous discretely. Undergraduate Texts in Mathe- 
matics. Springer, New York, 2007. 

[7] Christian Bey, Martin Hcnk, and Jorg M. Wills. Notes on the roots of Ehrhart polynomials. Discrete 
Comput. Geom., 38(l):81-98, 2007. 

22 



[8] Norman Biggs. Algebraic graph theory. Cambridge Mathematical Library. Cambridge University Press, 

Cambridge, second edition, 1993. 
[9] Benjamin Braun. An Ehrhart series formula for reflexive polytopes. Electron. J. Combin., 13(l):Note 

15, 5 pp. (electronic), 2006. 
[10] Benjamin Braun and Matthias Beck. Euler-Mahonian permutation statistics via polyhedral geometry. 

Preprint. 

[11] Katie Bright and Carla Savage. The geometry of lecture hall partitions and quadratic permutation 

statistics, to appear in FPSAC proceedings, 2010. 
[12] Sylvie Corteel, Carla D. Savage, and Herbert S. Wilf. A note on partitions and compositions defined by 

inequalities. Integers, 5(1):A24, 11 pp. (electronic), 2005. 
[13] J. William Davis, Erwin D'Souza, Sunyoung Lcc, and Carla D. Savage. Enumeration of integer solutions 

to linear inequalities defined by digraphs. In Integer points in polyhedra — geometry, number theory, 

representation theory, algebra, optimization, statistics, volume 452 of Contemp. Math., pages 79-91. 

Amer. Math. Soc, Providence, RI, 2008. 
[14] Eugene Ehrhart. Sur les polyedres rationnels homothetiques a n dimensions. C. R. Acad. Sci. Paris, 

254:616-618, 1962. 

[15] A. Elashvili, M. Jibladzc, and D. Pataraia. Combinatorics of necklaces and "Hermite reciprocity". J. 

Algebraic Combm., 10(2):173-188, 1999. 
[16] Christian Haase and Ilarion V. Melnikov. The reflexive dimension of a lattice polytope. Ann. Comb., 

10(2):211-217, 2006. 

[17] Takayuki Hibi. Dual polytopes of rational convex polytopes. Combinatorica, 12(2):237-240, 1992. 
[18] Matthias Koppe. Latte macchiato, version 1.2-mk-0.6. Available from URL 

http://www.math.ucdavis.edu/%7emkoeppe/latte/, 2009. 
[19] Jeffrey C. Lagarias and Giinter M. Ziegler. Bounds for lattice polytopes containing a fixed number of 

interior points in a sublattice. Canad. J. Math., 43(5):1022 -1035, 1991. 
[20] Mircea Mustata and Sam Payne. Ehrhart polynomials and stringy Betti numbers. Math. Ann., 

333(4):787-795, 2005. 

[21] Igor Pak. Partition identities and geometric bijections. Proc. Amer. Math. Soc, 132(12):3457-3462 
(electronic), 2004. 

[22] Sam Payne. Ehrhart series and lattice triangulations. Discrete Comput. Geom., 40(3):365-376, 2008. 
[23] Richard P. Stanley. Decompositions of rational convex polytopes. Ann. Discrete Math., 6:333-342, 1980. 

Combinatorial mathematics, optimal designs and their applications (Proc. Sympos. Combin. Math, and 

Optimal Design, Colorado State Univ., Fort Collins, Colo., 1978). 
[24] Richard P. Stanley. Enumerative combinatorics. Vol. 1, volume 49 of Cambridge Studies in Advanced 

Mathematics. Cambridge University Press, Cambridge, 1997. With a foreword by Gian-Carlo Rota, 

Corrected reprint of the 1986 original. 
[25] Giinter M. Ziegler. Lectures on polytopes, volume 152 of Graduate Texts in Mathematics. Springer- Verlag, 

New York, 1995. 



23 



Department of Mathematics, University of Kentucky, Lexington, KY 40506-0027 
E-mail address: benjamin.braun@uky.edu 

E-mail address: davis.robert@uky.edu 

E-mail address: ashley.liarrison0525@uky.edu 

E-mail address: jessica.inckiin@uky.edu 

E-mail address: jenna.noll@uky.edu 

E-mail address: clifford.taylor@uky.edu 



24 



